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fNl ■ canonical method based on Bogoliubov transformation is applied to calculate the pair 

. creation probability and the density number of created particles. Some particular 

cosmological models such as radiation dominated universe and Milne universe are 
discussed. For both cases the vacuum to vacuum transition probability is calculated 
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INTRODUCTION 



As we know there is no well defined theory of quantum gravity and the quantization 
of gravitational fields is one of unsolved problems in theoretical physics. Therefore the 
straightforward way to treat gravitational effects is the quantum field theory in curved 
space-time which proved most fruitful in describing interaction of matter with gravity and 
finding physically meaningful quantities [irfll. The quantum field theory in a curved space 
is the first approximation of quantum gravity where the gravitational field described by 
the metric of the space-time is treated as classical field and the matter is described by 
quantum fields which propagate in such a gravitational field. It is well known also that the 
most significant prediction of this theory is the phenomenon of particle creation from the 
vacuum. Generally in curved space-time it is difficult to define a physical vacuum state for 
the quantum field and when this vacuum state is defined in the remote past it is usually 
unstable so that it may differ from the vacuum state in the remote future, and spontaneous 
creation of particles occurs . 

The importance of particle production by gravitational fields comes from the diversity of 
its applications in contemporary cosmology; Particle creation effects could have consequences 
for early universe cosmology and could play a role in the exit from inflationary universe and 



in the cosmic evolution [6HlO|. The application of this phenomenon in black hole physics, 
such as Hawking radiation, is also well known [Tlj ] . 

In order to study the phenomenon of particle creation in gravitational fields we found in 
literature several techniques such as the adiabatic approach 
nalization method [l5-2(|, the Green function approach EH, 



technique [23l l24j. the "in" and "out" states formalism [25[ 



12j-|l4| , the Hamiltonian diago- 
, the Feynman path integral 
as well as the semiclassical 
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WKB approximation |29|-|33|. In addition we cite the technique based on vacuum to vacuum 



transition amplitude and Schwinger-like effective action [341 . 135 . 

The aim of this paper is to study the creation of scalar particles in a flat Robertson- Walker 
space-time provided with metric of the form 

ds 2 = dt 2 - a 2 (t) (dx 2 + dy 2 + dz 2 ) (1) 

In terms of conformal time 77 = J dt/a(t) this metric reads 

ds 2 = C (77) [drf - dx 2 - dy 2 - dz 2 ] (2) 

where the new scale factor C (77) is defined from a (t) as follows C (77) = a 2 (77) = a 2 [t (77)] . 
We choose for the scale factor the form 

C (77) = a + b tanh (A77) + c tanh 2 (A77) (3) 

where a, b and c are positive parameters. 

We can see that this form is the generalization of various particular cases found in liter- 
ature; When c = 0, we have a cosmological model with a 2 (77) = a + b tanh (A77) which has 
been widely studied (36 | -[38 | . With a particular choice of parameters a, b and c we get some 

models discussed in [39|, |40|. In addition, this universe becomes a radiation dominated one 

4 

when a = b = 0, c = and A — > 0. We can also make connection with a Milne universe 

2 

(i.e. a (t) = ait.) when c = 0, X — a±, b — a — ^ by making the change 77 — > 77 + ^and 
taking the limit e — > 0. 

In the first stage we introduce a scalar field propagating in Robertson- Walker space-time 
and we give two sets of exact solutions for the Klein Gordon equation. Next we use the 
relation between these two sets to determine the probability of pair creation, the density 
number of created particles and the vacuum persistence. Finally we discuss some particular 
examples. 

II. SCALAR FIELD AND KLEIN GORDON EQUATION 

Let us consider a scalar matter field $ with mass m subjected to the gravitational field 
described by the metric g^. The dynamics of this system is in general governed by the 
action 

S = J d 4 x [yS/ (-g^d^d;® - m 2 $$* - V (g) $$*)] (4) 

where V (g) is a function of invariant combinations of the metric tensor g^ u and its partial 
derivatives 42 

V(g)=bR + &R 2 GB , (5) 

where R is the Ricci scalar, £1 and £2 are two numerical factors and Rqb is a scalar which 
explains the Gauss-Bonnet coupling 

Rgb = R 2 — ^R^vR^ + R^vpaR IJ,upcr ■ (6) 
In such a case the Klein Gordon equation can be written in the form 
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,d„ (g^y/^d^) + (m 2 + V (g)) $ = 0. (7) 



In the present work we interest only on conformally coupled scalar field (i.e. £1 = ~ and 
£2 = 0). If we choose to work with conformal time rj which is convenient to the present 
coupling and we introduce a new field ip (x) so that 

$ W = /777 — (x) = ~J^=X (x) V? (v) , (8) 

where \ {%) has, in the case of flat space-time, the form of a plane wave \ {%) ~ ex P \ik.rj , 
we can obtain the simplified equation 



d 2 , 



2 (V) 



<p(ri) = 0, (9) 



with 

cu 2 (ri) = k 2 + m 2 C(r]). (10) 

Since equation ()9]) is of second order there are only two independent solutions and all other 
solutions can be expressed in terms of these two independent ones. Here we want to find 
two sets of independent solutions so that the two functions tpf n of the first set behave like 
positive and negative energy states at r\ — > — oo and the two functions (pf ut of the second 
set behave like positive and negative energy states at 77 — > +00. Therefore, before to look 
for exact solutions for equation let us first, examine their behavior when 77 — > ±00. We 
easily obtain 

^W=exp(-«eu m f/) (11) 
Vint (v) = exp {-ieu out r]) , (12) 

where e indicates the positive or the the negative frequency mode ( e = ±1) and u out and 
u in are given by 



oJin = -\A 2 + m 2 (a + c - b) (13) 
ooout = ^k 2 + m 2 (a + c + b). (14) 

Note that behaviors in f lTTj) and f|T2|) are the same as semi-classical solutions obtained from 
Hamilton- Jacobi equation because the space-time is asymtotically Minkowskian. 
Now in order to solve equation (Q we make the change 77 — > £, where 

f _ 1 + tanh ( X V) ( 15 ) 
The resulting equation that takes the form 
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is a Riemann type equation 43 



aci — a[ 1 — «3 — a' 3 \ d 



aidi . 
— a 2 a 2 + 



where 



1-^ 
1 



^(0 = o 



(17) 



«2 



—a 
—a 

1 - 



«2 



2A 

2A 
1 
2 



(18) 



A 2 



1 

4' 



with the condition ai + a[ + a 2 + a' 2 + a 3 + a' 3 = 1. 



Following [43| we can find for equation ([9]) several sets of solutions that can be written 
in terms of hypergeometric functions. Taking into account the behavior of positive and 
negative energy states we can classify our to sets as follows; for the "in" states we have 



and 



with 



and 



^(0 = ^r-a-o^ 

V AUJ in 

-c- 1 -W- . c- , .Win i 

V 2 + 1 ~ + ' 2 + % ~ ~ ! X ; ^ ' ( } 



'fin (0 = -rx=Z 2X (! - 2A 



w ± = o ( 21 ) 



(22) 



1 Am?c 

The factors (2a;j n ) -1 ' 2 and (2a; ou t)~ 1 ^ 2 are determined by the use of the following normaliza- 
tion condition 

V*k<Pk - VkVl = 2i, (23) 

which explains the conservation of the Klein Gordon particle current. 
For the "out" states we have 



V ^W out 



\ + i^ + i5,\ + i^-i5^ + i?f-,l-Z) (24) 
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and 



Pout (0 



2A 



out 
1 W_ 

2 _, T 



' 2 ~ - id; 1 - z— -; 



1-^ 



(25) 



Having shown how to solve Klein Gordon equation for spinless particle conformally cou- 
pled to gravitational field described by the metric (J2J), let us investigate the particle creation 
process. In the next section we quantify the scalar matter field and we use the solutions of 
the field equation to analyze the problem in question. 

III. FIELD QUANTIZATION AND PARTICLE CREATION 



In the first stage we write the field operator in it's Fourier decomposition 



^ ^ = 71 ' d k ^ ^ Xk ^ ak + Vk ^ x * k ^ 6 * 



(26) 



where, in canonical quantization formalism, the operators dk and bk verify the following 
commutation relation 



8 k-k' 



(27) 



With the help of the normalization condition ( 123|) we can find without difficulties the fol- 
lowing expression of the Hamiltonian associated with the scalar field system 



H 



E k (v) ( a k al + bfb k ) + F* k (rj) b k a k + F k (rj) a+bj 



with 



E k (v) = \<Pk (v)\ 2 + ^1 (rj) \ip k (rj)f 

Fk (v) = <pI iv) + u 2 k (v) d (v) ■ 



(28) 



(29) 
(30) 



Here we remark that H is not diagonal at any time. For rj — > ±oo, however, Fk (rj) = and 
H becomes diagonal. In this situation we can define two vacuum states \0 in ) and \0 out ) . 
The state |0j„) is an initial quantum vacuum state in the remote past with respect to a static 
observer and \0 out ) is a final quantum vacuum state in the remote future with respect to the 
same observer. This gives some vacuum instability which leads to particle creation. 

Now, as we mentioned above, in order to determine the probability of pair creation 
and the density of created particles we use the solutions of the Klein Gordon equation by 
considering the canonical method based the relation between "in" and "out" states. In other 
wards Bogoliubov transformation connecting the "in" with the "out" states can be projected 
in Fock space to be converted into relation between the creation and annihilation operators; 
therefore, the probability of pair creation and the density of created particles will be given 
in terms of Bogoliubov coefficients. 
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So to get connection between "in" and "out" modes let us use the relation between 



hypergeometric functions [43 



^ / r (w) r (w — v — u) _ , _ _ ^ x 

-F (tt, f ; w; ^) = ! \ >-F (u, v; u + v - w + 1; 1 - 

1 [w — u) 1 [w — v) 

fl _ r (7) r (u + v - w) 
1 * J i>)r» 

F (it; — to — u> — v — u + 1; 1 — £) (31) 

and the property 

F (it, v- w; = (1 - F (w — u,w — v;w; f) (32) 

to obtain the so-called Bogoliubov transformation connecting "in" and "out" states 

V>tn = a ¥tut + PVout ( 33 ) 
fin = P*<Ptut + a *fout, (34) 

where the Bogoliubov coefficients a and (3 are given by 

a= r(i-^)r(-^) (35) 



v " - 'IT - <-\/ — - 1 J 1 \~ - '— <-\/ — 

and 

fi = ,/^- r(i-iy)r(« T ) _ (M) 

V <*» r ^ + ; - _ ii^-i) r (i + ^ + - 1 

with \a\ 2 — |/3 1 2 = l.The relation between the creation and annihilation operators is then 



a out = a a in + f3*bf n (37) 
bt ut = P a in + a*b+. (38) 

For the process of particle creation the probability amplitude that we want to calculate is 
defined by 

A=(0 out l&outUout I in ) . (39) 

Taking into account that 

1 (3* 

bout = —kn + —at ut (40) 



we obtain 



f3* 

■A = (0 ou t \<loutbout \ Oin) — — (Oout I Oin) ■ (41) 

a* 
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The probability to create one pair of particles from vacuum is then 

(3* 2 



or 



Pcreat. \k) 

Using the following properties of the Gamma functions [43 

r(z+ 1) = zT(z 



and 



we arrive at 



\r(ix) 



T(- + ix 



TV 



x sinh tix 



71 



creat. 



(k) 



COsh 7TX 

cosh (2vr^) + cosh (2tt5) 
cosh (2tt^) +cosh(27r5)' 



(42) 

(43) 
(44) 

(45) 
(46) 



Let Ck to be the probability to have no pair creation in the state k.The quantity Ck {P cr eat. 
is then the probability to have only n pairs in the state k. We have 



Ck ~\~ CkPcreat. Ck (-Pcreat.) Ck (Pcreat.) ~\~ ■■■ 1 



(47) 



or simply 



Being aware of I 1 2 + 1-^1 



Cfc = 1 — Pcreat.- (48) 

1, we can find the vacuum persistence which reads 
i 2 



Ci 



a* 



cosh (2tt^) - cosh (2tt^) 



cosh (2tt^) + cosh (2tt5) 
Another important result is the density number of created particles 

n (k) = (0 in {agouti in ) = \(3\ 2 

A simple calculation gives 



(49) 



n (k) 



cosh {2tt^A + cosh 



_ , 4m 2 c 1 
TT\I — 1 



cosh (2tt^) - cosh (2vr^) 
Here we note that the density number of created particles can be written as 

n (k) 





2 


a 


- 1 








(50) 



(51) 
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and for large frequencies n (k) becomes a thermal Bose-Einstein distribution 

n(k) = * , - ■ (52) 

exp {—LJin) ~ 1 

From equation ( |50|) one can see that the pair creation process is more important for low 
impulsion and n (k) — > when k — > oo (i.e. k » 1). Also, for a — b + c ^ and a + b > c, 
it is clear that n (k) — > when m >> 1. 



IV. SPECIAL CASES 

In order to illustrate our calculation let us in this section discuss some particular cos- 
mological models such as radiation dominated universe and the Milne universe which are of 
interest and have been studied by other authors and by using other techniques. We start by 
radiation dominated universe. 



A. Radiation dominated universe 

4 

Considering the case when a = b = and c = an d taking the limit A — > 0, we make 
connection with the well-known phase of radiation dominated universe 



and 



a (t) = a Q Vt. 

In this case the probability to create a pair of particles will be 



(53) 
(54) 



exp 7T 



V rad . (k) = lim ■ 



4m 2 c 
A 2 



1 



2-kuj 
A 



1 + exp ( 7i 



im?c i 2-7TOJ 
1? 1 ~ 



Taking into account that 



(55) 



'4m 2 c 
n\l — 1 



2ttu ai I AX 2 k 2 



A 



1 a 2 



m 2 ag 2 A 



X 4 



man 



(56) 



we get 



V r ad. (k) 



exp 




2k 2 ' 
mai 




1 + exp 


2k 2 
mag 



The density number of created particles will be 



(57) 



n (k) = exp 



mat 



(58) 
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This results coincide exactly with those found in literature 23 . 
For the vacuum to vacuum probability we can write 



V vac ^ vac = exp (-2 Im S eff ) 
= \{C k 

k 

= IJ ex P t~ ln (! + (T )] 

k 

= exp — In (1 + a) 

k 

where a = exp ^— ■ Consequently, we have 

2ImS e// = £>(l + <r). 



Expanding the quantity ln(l + a), we get 

d 3 k 



2Im£ 



eff 



(2*r ^ 



v (-l) r 



,n+l 



exp — 717T- 



2A; 2 



By doing integration over momentum k 



2Im£ 



E— 



,n+l 



n=l 



5 

n 2 



and summing over n, we obtain the following result 



2ImS e// = 0.867 



mat 



(59) 



(60) 



(61) 



(62) 



(63) 



It is clear that the effect becomes important as soon as the factor ao approaches the critical 
value 



a 2 ~ 

cr 



and the pair production rate is of order of 



m 

3 



B. Milne universe 

Now we consider the example of Milne universe defined by the scale factor a (t) = ait. In 
term of conformal time this scale factor reads 



C (rf) = a 2 (77) = a\ exp (2a\rj) 



(65) 
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As is mentioned above this case may be obtained from ([3]) by considering the case c = 0, 

a 2 1 

X = a\, b = a = making the change 77 — )• 77 + -|^and taking the limit e — > 0. In these 
conditions the probability to have one pair creation will be 



7W=exp[-^ 
ai 



(66) 



which is in agreement with the same probability calculated by other authors [23 
The vacuum persistence is 

( 2tt, 

Cfc = 1 — exp k 

V a i 

and the vacuum to vacuum probability takes the form 



(67) 



exp (—2 Im S, 



eff) 



exp 



from which we draw the imaginary part of the effective action 

2Im% = 



Making the Taylor expansion of the logarithm function 

d 3 k 



21m S e f f 
and summing over n we obtain 



(27T) 



V-expf-n^ 



(68) 



(69) 



(70) 



21m S. 



1.2 x 



4vr 4 



Here we remark that the pair production rate is of order of a\. 



(71) 



C. Other particular cases 

Since the creation of superheavy particles with the mass of the Grand Unification scale 
in the early Universe has many important cosmological consequences 4J-|46| we want to 
discuss in this paragraph cases with m » 1. In general when m » 1, we have shown that 
n (k) — > 0. However, when a — b + c = or c > a + b we remark that n (k) — > C st > even 
if m » 1. This implies that it is possible to create superheavy particles in this model. 

Let us mention that if we put c = 0, we obtain the same results as those found in 37-4(3 
and when a = (^) ; b = ^f-; c = (^) our results coincide with those of 39, 4^ . 
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V. CONCLUSION 

In this paper we have studied the creation of scalar particles in some flat Robertson- 
Walker space-times. We have considered the canonical method based on Bogoliubov trans- 
formation connecting the "in" with the "out" states. We have given two sets of exact 
solutions for the Klein Gordon field equation and we have used these solutions to expressed 
the probability of pair creation and the density of created particles in some closed forms. 
Then we have discussed some particular cosmological models where our results become the 
same as those obtained by other authors. For the radiation dominated universe and the 
Milne universe we have calculated the vacuum to vacuum transition probability and we 
have extracted the nonvanishing imaginary term of the effective action that means that cre- 
ated particles are real and not virtual ones. We have shown also that creation of superheavy 
particles may be important in the case when a — b + c = or c > a + b. 
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